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Abstract. The high-frequency peaked BL Lac PKS 2155—304 at redshift z=0.117 has been detected with high 
significance (~45<r) at energies greater than 160 GeV, using the H.E.S.S. stereoscopic array of imaging air- 
Cherenkov telescopes in Namibia. A strong signal is found in each of the data sets corresponding to the dark periods 
of July and October, 2002, and June-September, 2003. The observed flux of VHE gamma rays shows variability 
on time scales of months, days, and hours. The monthly-averaged integral flux above 300 GeV varies between 10% 
and 60% of the flux observed from the Crab Nebula. Energy spectra are measured for these individual periods of 
data taking and are characterized by a steep power law with a time-averaged photon index of T = 3.32 ± 0.06. 
An improved \ 2 P er degree of freedom is found when either a power law with an exponential cutoff energy or a 
broken power law are fit to the time-averaged energy spectrum. However, the significance of the improvement is 
marginal (~2<r). The suggested presence of features in the energy spectrum may be intrinsic to the emission from 
the blazar, or an indication of absorption of TeV gamma rays by the extragalactic infrared background light. 
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1. Introduction 

The high-frequency peaked BL Lac (HBL) object PKS 
2155-304 (z = 0.117) was discovered by the HEAP 1 X- 
ray satellite ^Griffiths et al, 1979| |Schwartz et al, 1979) 1 
at a position consistent with a poorly localized detec- 
tion by the Ariel V satellite (Cooke et al, 1978) 
PKS 2155—304 is now known to be one of the 
brightest extragalactic X-ray sources in the sky 
and as a result other X-ray satellites, including 
RPSAT, BeppoSAX, RXTE, and Chandra (see, e.g., 



Brinkmann et al, 1994 


Chiappetti & Torroni, 1997 


Vcstrand & Srcekumar, 


1999 


|Nicastro et al., 2002J, 



have detected it on a regular basis. PKS 2155—304 is well 
studied and has a history of strong broad-band variabil- 
ity. It is associated with a compact, flat-spectrum radio 
source, and exhibits an essentially featureless continuum 
from radio to X-ray frequencies. The maximum power 
emitted by PKS 2155-304 is between the UV and soft 
X-ray range, and it is the brightest BL Lac detected 
in the UV regime (Wandel & Urry, 1991) Gamma-ray 



emission in the energy range 30 MeV to 10 GeV was 
detected from PKS 2155-304 by the EGRET detector 
aboard the Compton Gamma Ray Gbservatory satellite 



(Vestrand et al, 1995)| The EGRET observations indi- 
cated a hard energy spectrum with integral power-law 
photon index of 1.71 ± 0.24. All of these characteristics 
make PKS 2155-304 a likely emitter of >100 GeV gamma 
rays. In fact, it was detected at energies greater than 
300 GeV in 1996 and 1997 by the University of Durham 



Mark 6 Telescope as reported in Chadwick et al. (1999a) 
However, the object was not detected by the same instru- 
ment in 1998 



(Chadwick et al, 1999b) 



Furthermore, 

PKS 2155-304 was not detected by the CANGARPP 
experiment during observations made in 1997 
(Roberts et al, 1999)1 1999 



(Nishijima et al, 2001) 



2000 or 2001 |(Nakase et al, 2003)| The reported upper 
limit from CANGARPP in 1997 is consistent with the 
flux detected by Durham, and the lack of a detection of 
PKS 2155—304 at TeV energies in subsequent years is 
not inconsistent with the fact that such emission from 
AGN is known to be highly variable. 

The present confirmation of the detection of VHE 
gamma-rays from PKS 2155—304 by the High Energy 
Stereoscopic System (H.E.S.S.) along with observations at 
other wavelengths will yield considerable insight into the 
mechanisms for TeV gamma-ray emission from blazars, 
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help delineate the differences between low-frequency 
peaked BL Lacs and HBLs, and will assist in the un- 
derstanding of the extragalactic infrared background light 
(EBL). This latter point is strengthened by the fact that 
PKS 2155—304 is the second most distant object de- 
tected at TeV energies and therefore its energy spec- 
trum may exhibit characteristics, such as steepening of 
the spectrum and a cutoff, caused by absorption of VHE 
gamma rays via pair production on the EBL. The most 
distant object detected at TeV energies, 1ES 1426+428 
(z = 0.129), showed features in its energy spectrum which 
have been interpreted as consequences of such absorption 



(Aharonian et al, 2002a) In addition, variability studies 
using detailed light curves obtained from observations of 
PKS 2155—304 can place strong constraints on the phys- 
ical modeling of blazars. 

2. H.E.S.S. Detector 

The H.E.S.S. experiment has been operating since June, 
2002, in the Khomas Highlands of Namibia (23° 16' 18" 
S, 16° 30' 1" E, 1835 m above sea level). The detec- 
tor consists of a system of four imaging air-Cherenkov 
telescopes in a square of 120 m side. Each individ- 
ual telescope is an alt-az mount Davies- Cotton reflector 



(Davies & Cotton, 1957) with a flat-to-flat width of 13 m, 
and has a camera mounted at the focal length of 15 m. 
The total mirror area of each telescope is 107 m 2 , seg- 
mented into 382 individual round (60 cm diameter) front- 
aluminized glass mirrors. The H.E.S.S. cameras provide 
a 5° field of view and contain 960 individual photomulti- 
plier (PMT) pixels subtending 0.16° each, with Winston 
cone light concentrators. The camera is modular in de- 
sign, housing 60 drawers of 16 PMTs each, and contains 
all the necessary electronics for operation, triggering, and 
readout. The trigger electronics divide the camera into 
overlapping 64 PMT sectors (4 adjacent drawers) with a 
trigger requirement that a sector has a minimum number 
of pixels with a signal above a threshold in photoelectrons 
(PEs) coincident in an effective ~1.3 ns trigger window. 
Upon receiving a camera trigger, the signal stored in ana- 
log memories from each of the PMTs, sampled at 1 GHz, 
is integrated within a 16 ns window. Once a camera has 
triggered, a signal is sent out via an optical fiber to a cen- 
tral trigger system (Funk et al, 2004) which allows for 
a multiple telescope coincidence requirement. More de- 



tails on H.E.S.S. can be found in Bernlohr et al. (2003)| 
Hofmann (2003)| |Vincent et al. (2003) ~ 



3. Observations 

The H.E.S.S. observations of PKS 2155-304 were made 
while the system was under construction, therefore the 
data cover various telescope configurations and trigger 
criteria. Observations in 2002 were made with just one 
telescope. In January 2003, a second telescope was added 
to the array and observations in 2003 consist mostly of 
a two-telescope configuration. The exception to this is in 
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Table 1. Shown are the dark periods in which PKS 2155—304 was observed, the configuration of H.E.S.S. during those 
observations, the telescope multiplicity required to trigger the array, the number of telescopes available, the sector 
trigger requirement (number of pixels above a threshold in PEs), the observation mode, the number of on-source runs, 
the dead time corrected observation time, the system rate, the dead time percentage, and the mean zenith angle of 
the observations (Z b s ). 



Dark 




Ntel 


Sector 


Obs. 


No. 


Obs. Time 


Syst. Rate 


Dead Time 


Zobs 


Period 


Configuration 


Mult., Tot. 


Trigger 


Mode 


Runs 


[hrs live] 


[Hz] 


% 


n 


07/2002 


Mono (30 Drawers) 


1,1 


4 > 6.7 


On - Off 


8 


2.5 


130 


24 


15 


10/2002 


Mono (52 Drawers) 


1,1 


4 > 6.7 


On - Off 


16 


4.3 


230 


37 


16 


11/2002 


Mono 


1,1 


4 > 6.7 


On - Off 


4 


0.8 


250 


53 


35 


06/2003 


Software Stereo 


2,2 


4 > 6.7 


Wobble 


32 


10.8 


60 


23 


13 


07/2003 


Hardware Stereo 


2,2 


3 > 5.3 


Wobble 


55 


22.1 


100 


7 


24 


08/2003 


Hardware Stereo 


2,2 


3 > 5.3 


Wobble 


45 


19.6 


100 


7 


14 


09/2003 


Hardware Stereo 


2,2 


3 > 5.3 


Wobble 


2 


0.9 


120 


8 


8 


09/2003 


Hardware Stereo 


2,3 


3 > 5.3 


Wobble 


6 


2.1 


170 


12 


18 



September 2003, when a third telescope was added to the 
array. The fourth and final telescope was added to the 
array in December 2003, after the end of the observation 
season for PKS 2155-304. 

The trigger for H.E.S.S. has similarly evolved. During 
the observations in July, October, November 2002, and 
June 2003, the sector trigger consisted of a requirement 
for 4 pixels each with more than 6.7 PEs. However, the 
July and October 2002 observations were made with only 
the inner 30 and 52 drawers, respectively, participating in 
the trigger while the entire camera was used for image col- 
lection. For all later observations the cameras were used 
fully. The system level trigger was installed in July 2003. 
Before this time two-telescope data were taken with each 
telescope separately, and the stereo multiplicity require- 
ment (two telescopes) was performed off-line ("Software 
Stereo") using GPS time stamps. After the installation 
of the central trigger system the stereo multiplicity re- 
quirement was performed in the hardware ("Hardware 
Stereo"). As the system rate was considerably lower due to 
a two-telescope multiplicity requirement, the sector trig- 
ger requirement could be lowered to 3 pixels with more 
than 5.3 PEs, allowing for a lower energy threshold. This 
increased the individual telescope rates, from ^250 Hz to 
^800 Hz, while maintaining a reduced system dead time. 

Table Q] gives details of the observations of PKS 
2155—304 by H.E.S.S. which pass conservative run selec- 
tion criteria. These run selection criteria remove runs for 
which the sky was not clear and where the telescopes were 
not operating within specified requirements (e.g., trigger 
rate stability, more than 95% of the camera pixels op- 
erational). All the single telescope observation runs were 
taken in On- Off mode. Here, data were taken on source 
for 25 minutes, preceded or followed by off-source observa- 
tions, used to estimate the background. These off-source 
observations differ by ±30 minutes in right ascension (RA) 
from the actual source position. Given an approximately 5 
minute (reduced to 2 minutes after 2002) transition time 
for the data acquisition system, this ensures that the on 
and off runs have identical azimuth and altitude profiles. 



The stereo data were taken in 28 minute runs using Wobble 
mode. In this mode, the source direction is positioned 
±0.5° in declination relative to the center of the field of 
view of the camera during observations. The sign of the 
offset is alternated in successive scans to reduce system- 
atic effects. Due to the large field of view of the H.E.S.S. 
cameras, use of Wobble mode allows for both on-source ob- 
servations and simultaneous estimation of the background 
induced by charged cosmic rays, since the background can 
be estimated from different regions in the same field of 
view. This eliminates the need for off-source observations 
and therefore doubles the amount of time available for 
on-source observations. 



4. Technique 

4.1. Event Reconstruction 

The analysis of the data passing the run selection crite- 
ria proceeds in the following steps: First the images are 
calibrated (Aharonian et al., 2004) and then "cleaned" 
to remove noise from the image. The image cleaning is 
done using a two-stage tail-cut procedure which requires 
a pixel to have a signal greater than 10 PEs and a neigh- 
boring pixel to have a signal larger than 5 PEs. Also pixels 
greater than 5 PEs are included if they have a neighbor 
greater than 10 PEs. After this image cleaning is per- 
formed, the moments of the shower image are parame- 
terized using a Hillas-type analysis (Hillas, 1985) The 
shower geometry is reconstructed using stereoscopic tech- 
niques (Aharonian et al., 1999b) with a typical angular 

TO m in 



resolution of ^0.1° and an average accuracy of 
the determination of the shower core location. To ensure 
that the analyzed images are not truncated by the edge of 
the camera, only images which pass a distance cut requir- 
ing the image center of gravity to be less than 2° from the 
center of the camera are used in the reconstruction. For 
the stereo configuration a minimum of 2 telescopes meet- 
ing this criterion are required for a successful reconstruc- 
tion. In addition, at least two telescopes are each required 
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Table 2. The selection cuts applied to the data and the percentage of 7-ray and background events retained by those 
cuts. 







MRSL 


MRSL 


MRSW 


MRSW 


C.O.G. 


Size 


Length/Size 


Ntel 


2 


7 


BG 


Configuration 


Ntel 


min 


max 


min 


max 


max 


min 


max 


min 


max 


% 


% 






M 


M 


M 


w 


[mrad] 


[PE] 


[mrad/PE] 




[deg 2 ] 






Mono 


1 


-5.0 


3.4 


-4.2 


1.0 


35 


None 


0.016 


1 


0.035 


20 


0.012 


Stereo(Soft.) 


2 


-2.2 


3.2 


-10.0 


1.6 


35 


75 


None 


2 


0.030 


72 


0.036 


Stereo (Hard.) 


2 


-1.8 


2.0 


-10.0 


1.0 


35 


55 


None 


2 


0.025 


48 


0.029 


Stereo (Hard.) 


3 


-1.8 


2.0 


-10.0 


1.0 


35 


55 


None 


2 


0.025 


54 


0.041 



to exceed a minimum total signal to ensure that the images 
are well reconstructed. The analysis techniques are similar 
for the mono-configurations, although a necessary differ- 
ence arises in the determination of the shower geometry. 
In the mono-telescope case, the impact parameter and di- 
rection of the shower axis are estimated directly from the 
image. To do this, the angular distance, along the major 
axis of the reconstructed ellipse, from the image center 
of gravity to the point of origin of the event is estimated 
using the image parameter construction length/log(size) 



(Masterson, 2000) 



4.2. Background Rejection 

After the event reconstruction, the much more numerous 
cosmic-ray background events are rejected using cuts on 
mean reduced scaled width (MRSW) and length (MRSL) 
parameters. These parameters are defined as the mean of 
the difference in standard deviations for each telescope of 
the width (length) observed in the image from that which 
is expected from gamma-ray simulations (<width> and a) 
based on image intensity, reconstructed impact parameter 
and zenith angle of observations. The equation for MRSW 
is as follows: 



MRSW = ^ E^o 



<width> i 



(1) 



and is similar for MRSL. Images from muons are rejected 
by requiring events be coincident in at least 2 telescopes in 
the stereo configuration or by using a cut on length divided 
by size of the image, after image cleaning, for the mono- 
configuration (Catanese et al, 1995) A cut on 2 , the 



square of the angular difference between the reconstructed 
shower position and the source position, is applied and is 
equivalent to placing the data into a round bin centered on 
the source position. All the cuts (shown in Table0) are op- 
timized a priori (simultaneously) to yield the maximum 
expected significance per hour of observation. The data 
used for optimization consists of Monte Carlo gamma-ray 
simulations at a zenith angle of 20° with a Crab-like en- 
ergy spectrum (spectral index=2.6, 1/2 Crab flux) and 
real off-source data. The percentage of gamma-rays and 
cosmic rays retained by the cuts is also shown in Table El 
The significance expected (and observed) is not strongly 
dependent on the exact values of the cuts used. For each 
of the various hardware configurations a different set of 



cuts is used, all optimized on Monte Carlo simulations 
and comparable off-source data sets. 

4.3. Background Estimation 

For mono-telescope observations the background is esti- 
mated from the number of events passing the selection 
cuts at the center of the field of view during off-source 
observations. The background is then normalized using 
the ratio of dead time corrected observation time. For the 
stereo ( Wobble) observations the background is estimated 
using all events passing cuts in a ring around the source 
location. The distance from the source position to the ra- 
dial center of the ring is 0.5°, and the width of the ring is 
adjusted such that the area of the ring is approximately 
seven times the area of the on-source region. The normal- 
ization of the off-source event total is corrected for the 
radial acceptance of the camera. The use of a larger back- 
ground region in this ring technique reduces the relative 
statistical error on the background measurement. After 
the number of on-source and background events passing 
the selection cuts are determined, the significance of the 
excess is calculated following the method of Equation (17) 
in|Li fc Ma (1983) 



5. Detection of PKS 2155-304 

PKS 2155—304 is strongly detected in all the dark periods 
in which it was observed with the exception of November 
2002 where the exposure was less than one hour. Table |21 
shows the results of the H.E.S.S. observations for each of 
the individual dark periods. The total significance of the 
excess for all observations is 44.9<r. Figure fj] shows the on- 
source and normalized off-source distributions of 9 2 for all 
observations in the two-telescope hardware stereo config- 
uration 1 . The background is flat in 9 2 as expected, and 
there is a clear excess at small values of 6 2 corresponding 
to the observed signal. Figure |U shows a two-dimensional 
sky map of the excess observed in the direction of PKS 
2155—304 for the same configuration. The bins are not cor- 
related and represent the actual distribution of observed 
gamma-rays (and some background events) on the sky. A 



1 This configuration represents the major share of the data; 
the shapes of the respective curves are similar for the other 
configurations. 
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Table 3. Shown are the number of on-source and off- 
source events passing the cuts, the normalization for 
the off-source events, the observed excess from PKS 
2155—304, and the significance of the excess for each of 
the dark periods. 



Dark 
Period 


Ntels 


On 


Off 


Norm 


Excess 


Sig. 

M 


07/2002 


1 


637 


234 


1.063 


388 


13.0 


10/2002 


1 


865 


519 


1.065 


312 


8.2 


11/2002 


1 


80 


90 


0.896 


-0.65 


-0.1 


06/2003 


2 


1396 


4619 


0.152 


694 


21.1 


07/2003 


2 


3169 


13643 


0.147 


1164 


22.1 


08/2003 


2 


3369 


12955 


0.147 


1459 


27.7 


09/2003 


2 


185 


697 


0.148 


82 


6.7 


09/2003 


3 


1005 


3946 


0.147 


425 


14.7 


Total 
Total 


Stereo 
All 










43.8 
44.9 



•a - 29 

0) 

o 

£-29.5 



-30 



-30.5 



-31 



-31.5 



22.05 22.00 21.95 21.90 

RA (hours) 





160 




140 




120 




100 




80 




60 




40 




20 









-20 




0.25 0.3 
e 2 (degree 21 



Fig. 1. The distribution of 6 2 for on-source events (points) 
and normalized off-source events (shaded) from observa- 
tions of PKS 2155—304 in the two-telescope hardware 
stereo configuration. For this plot, the background is taken 
from a region of comparable area located at the opposite 
offset position in the field of view. This causes the signif- 
icance shown to be reduced due to the smaller area used 
for the background estimation. 

fit of the peak in Figure [21 to a two-dimensional Gaussian 
finds the shape of the excess to be characteristic of a point 
source. The fit peak is located (Ara = 9 ± 12 sta t arc- 
sec, Adec = 14 ± 10 s tat arcsec) from the position of PKS 
2155—304, consistent with the position of the blazar as 
expected (systematic errors on the pointing are less than 
20 arcsec in RA and declination (Gillessen, 2003) I. 



Figure H shows the distributions of MRSW, MRSL, 
and 9 2 for the observed excess after all cuts are applied 
(excluding both the MRSW and MRSL cuts in (a) and (b), 
and the 9 2 cut in (c)), as well as the corresponding distri- 
butions for Monte Carlo gamma-ray simulations, for the 
entire two-telescope hardware stereo data set. The agree- 



Fig. 2. A two dimensional distribution of excess events 
observed in the direction of PKS 2155—304 for the two- 
telescope hardware stereo configuration. 

ment is good and shows that the simulations accurately 
characterize the detector response and that the detector 
performs as expected. In particular, the agreement of the 
9 2 distribution demonstrates that the point spread func- 
tion (width<0.1°) is well understood. 

6. Energy Spectrum 

6.1. Technique 

The energy of observed gamma rays, Eg t , is calculated us- 
ing the mean of the energies estimated for each telescope 
with a typical event resolution of ~15%. Each of the in- 
dividual energy estimates use the image size, impact pa- 
rameter 2 of the event, and the zenith angle of observation 
(Z) for each telescope. The energy estimates are based on 
results from Monte Carlo gamma-ray simulations. For the 
determination of the energy spectrum the data are loga- 
rithmically binned (8 bins per decade) in energy. However, 
only bins for which the average bias in Efi t is less than 
10% are used, effectively placing a lower energy threshold 
on suitable events. This is done to eliminate any system- 
atic effects that might arise because energy estimates for 
events near the trigger threshold yield, on average, too 
large an energy. The range of bins used varies for each in- 
dividual run due to the zenith angle of observation and re- 
sults in the bins at lower energies having differing amounts 
of exposure. For each event passing cuts the effective area, 
A c ff(Efit,Z), is individually determined. The influence of 
the finite energy resolution has been taken into account in 
the effective areas which are fit as a function of Z and Egt 



2 The distance of the center of gravity of the image from 
the center of the field of view is used in the mono-telescope 
configuration. 
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(a) (b) (c) 

Fig. 3. The distributions of a) mean reduced scaled width (MRSW), b) mean reduced scaled length (MRSL), and c) 
9 2 for the background subtracted signal after selection cuts (excluding both the MRSW and MRSL cuts in (a) and 
(b), and the 6 2 cut in (c)) from observations of PKS 2155—304 (points) and similarly for Monte Carlo gamma-ray 
simulations (spectral index = 3.35, zenith angle = 20°) in the two-telescope hardware stereo configuration (line). 



Table 4. The results of power law fits to the differential energy spectrum of PKS 2155—304 for each dark period. 
Shown are I , T, the \ 2 an d degrees of freedom (d.o.f) for the fit, the integral flux above 1 TeV (I(>1 TeV)), I(>300 
GeV), I(> E t h) where E t h is the post-selection cuts energy threshold of the observations, and the percentage of these 
fluxes relative to the respective values measured by H.E.S.S. from the Crab Nebula. Only the statistical errors are 
shown. 



Dark 
Period 



Io 



X 2 (d.o.f.) I(>1 TeV) 



m 2 s TeV 



Crab 

% 



I(>300 GeV) 



Crab 

% 



Eth 
[GeV] 



I(>Eth 

1-7 



Crab 

% 



07/2002 


15.6±2.1 


2.84±0.24 


7.0 (7) 


85.1±16.1 


43 


7.77±1.63 


56 


300 


7.77±1.63 


56 


10/2002 


6.36±1.75 


3.10±0.46 


0.5 (4) 


30.4±10.6 


15 


3.79±1.63 


28 


300 


3.79±1.63 


28 


06/2003 


2.42±0.28 


3.56±0.17 


4.0 (5) 


9.45±1.27 


4.8 


2.06±0.37 


15 


260 


2.97±0.59 


17 


07/2003 


1.75±0.18 


3.26±0.11 


3.9 (7) 


7.73±0.87 


3.9 


1.18±0.15 


8.6 


170 


4.27±0.74 


12 


08/2003 


1.84±0.18 


3.36±0.09 


7.4 (7) 


7.78±0.82 


3.9 


1.34±0.16 


9.7 


170 


5.13±0.80 


15 


09/2003 


2.40±0.41 


3.42±0.15 


2.9 (6) 


9.91±1.79 


5.0 


1.82±0.38 


13 


160 


8.35±2.27 


22 



for each of the various observation configurations similar 



to the method in Mohanty et al. (1998) For each bin, the 
sum of 1/ A c ff, is calculated for both on-source (Oi) and 
off-source {Bi) events using only events from the appro- 
priate runs. Each bin in the spectrum then contains: 



F, 



{Oi-Bi. 
AEiti 



(2) 



where and AE^ are the exposure time and width in 
energy respectively for that bin. As the collection area 
depends weakly on the spectral shape assumed, an it- 
erative procedure is followed, where the initial spectral 
shape is assumed to be flat in vF v (differential power- 
law index=2.0). When deviations from this hypothesis are 
found the collection area is recalculated using the previous 
result and the spectrum is determined again. This proce- 
dure continues until convergence (typically 1-2 iterations) . 

6.2. Dark Period Spectra 

Using the observed excess of gamma rays, a differential 
spectrum for each darkness period during which PKS 



2155—304 was detected is generated, and is shown in 
Figure These spectra are fit to a power law: 



dN_ _ j I E 
dE lo \ E„ 



(3) 



where I Q is the differential flux normalization, E is the en- 
ergy at which the flux is normalized (1 TeV), and T is the 
spectral index of the source. The dashed lines in Figure 0] 
show the best \ 2 fit to the data of Equation © for the 
corresponding spectrum. Table 0] shows the best fit re- 
sults, integral flux (above several thresholds) determined 
from the fit parameters, and comparisons of these inte- 
gral fluxes to those measured by H.E.S.S. from the Crab 
Nebula. Due to the steeper spectral index determined for 
PKS 2155-304 (r crab = 2.63 ± 0.05), the percentage of 
the respective Crab fluxes in an individual darkness pe- 
riod is larger for lower thresholds. Although the level of 
VHE emission from AGN is known to be variable, the 
integral flux above 300 GeV measured for each of the ob- 
servation period s shows reasonable agreem ent with the re- 
sult reported by |Chadwick et al. (1999a)l I(>300 GeV) = 
(4.2±0.75 s tat±2.0 syst )xlO~ Y m _z s _i . In addition to the 
method described above, an alternative forward-folding 
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Fig. 4. The differential energy spectrum measured from PKS 2155—304 for the dark periods of a) July 2002, b) 
October 2002, c) June 2003, d) July 2003, e) August 2003, f) September 2003. The dashed lines represent the best 
X 2 fit of the data to a power law. 



method, used by CAT (Piron et al, 2001) was also ap- 
plied to quantify the individual spectra and yields very 
similar results to those shown in Table 0] 



6.3. Time-Averaged Spectrum 

A time-averaged energy spectrum from all the stereo 
(2003) data is shown in Figure[S] The mono-telescope data 
set is not included for two reasons: The average flux level 
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Fig. 5. The time-averaged energy spectrum of PKS 
2155-304 for the 2003 (stereo) data set (top). The 
90% confidence level upper limits are determined by the 
method of Helene (1983) and are not used in the fits. The 
lines represent the best \ 2 fit to a power law (dot-dash), 
a power law with an exponential cutoff (dashed), and a 
broken power law (dotted). The bottom figure shows the 
residuals, (measurement — fit) divided by the statistical 
error, of the best fits to a power law (circles), power law 
with exponential cutoff (squares), and broken power law 
(triangles). 



in 2002 is considerably higher than the rest of the data set. 
More importantly, the energy is more poorly determined 
for the mono-telescope data set and thus the determina- 
tion of the spectrum is more susceptible to systematic er- 
rors. A fit of the stereo data set to a power law yields 
r=3.32±0.06, I Q = (1.96±0.12) x 10" 8 m" 2 s" 1 TeV" 1 , 
and a \ 2 of 10.8 for 7 degrees of freedom. Only the statis- 
tical errors are presented. Although the \ 2 f° r the power 
law fits to the spectra for individual dark periods is reason- 
able, the x 2 f° r the time-averaged spectrum corresponds 
to a probability of 0.148. In addition, the power law fit to 
the time-averaged spectrum is above the 90% upper limit 
on the differential flux at 3 TeV (not included in the fit). 
Since this information indicates that a power law is not a 
good characterization of the data, the energy spectrum is 



fit to two alternative hypotheses. The first is a power law 
with an exponential cutoff: 



dN_ _ t I E 
dE 2 ° [ E„ 



(4) 



where E cu t is the cutoff energy, and I , E , and V are 
as defined above. The fit yields an improved \ 2 of 6.2 
for 6 degrees of freedom (probability = 0.40), and the 
fit curve falls below all upper limits at higher energies. 
The best fit values are: r=2.90t£; 2 ,3, v 4 - ^!®) x 10~ 8 
m~ 2 s" 1 TeV -1 , and E cut = 1.4±£? TeV. The large errors 
(only statistical errors are shown) on the fit values are 
the result of strong correlations between the parameters 
used to characterize the spectral shape. Although the \ 2 
resulting from the fit is improved, it cannot be concluded 
that the fit is significantly better since an F-test shows 
there is a probability of only 0.92 that the power law with 
an exponential cutoff is a better characterization of the 
data. 

The second hypothesis is a broken power law, which 
follows Equation J2J for energies less than -Ebreak, and the 
following for energies greater than inbreak: 



dN _ j ( j 
dE — ±0 \~ 



(5) 



where -Ebreak is the energy of the break point in TeV, T\ 
and T2 are the spectral indices, and I and E are as de- 
fined above. The fit yields a \ 2 °f 5-1 for 5 degrees of 
freedom (probability = 0.40), and the fit curve does not 
fall below the upper limit at 3 TeV. The best fit values 
are: r 1 =3.15l°; 1 °, J = (2.4±g;f) x 10~ 8 m~ 2 s" 1 TeV" 1 , 

r 2 =3.79±o:27> and ^brcak = 0.7 ± 0.2 TeV. Only the sta- 
tistical errors are shown. An F-test shows there is a prob- 
ability of 0.85 that the broken power law is a better char- 
acterization of the data than the simple power law and 
again it cannot be concluded that the fit is significantly 
better. 

The possible presence of features, such as an exponen- 
tial cutoff or break, in the energy spectrum could be the 
result of absorption of TeV gamma rays by the extragalac- 
tic infrared background. However, it may also be intrinsic 
to the source spectrum, or the result of absorption of pho- 
tons close to the source, or both. The distinction can only 
be made using a larger sample of TeV-emitting AGN at 
different redshifts, or perhaps by detailed modelling of the 
broad-band spectrum using data gathered simultaneously 
at other wavelengths. 

7. Flux Variability 

7.1. Monthly Time Scale 

Figure [5] shows the integral flux greater than 300 GeV 
(top) and fit spectral index (bottom) of PKS 2155—304 
for each dark period. The dashed line in each plot repre- 
sents the results from a \ 2 fit of the data to a constant. 
The fit of I(>300 GeV) versus observation period yields 
a x 2 °f 24.0 for 5 degrees of freedom, corresponding to a 
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Fig. 6. The integral flux above 300 GeV and spectral in- 
dex versus dark period determined from H.E.S.S. obser- 
vations of PKS 2155—304. The dashed line represents the 
best x 2 fit of the data points to a constant. Only the sta- 
tistical errors are shown. 



2 x 10 x 2 probability. The data are clearly inconsistent 
with a constant flux. Although only the statistical errors 
are used in the fit, the conclusion is unchanged even if sys- 
tematic errors (~20% for the integral flux) are included. 
Since PKS 2155-304 was not detected in November 2002, 
a 99% upper limit of I(>300 GeV) = 3.0xlQ- 7 m ~ 2 s" 1 

assum- 



was determined by the method of Helene (1983) 
ing the 2003 time-averaged spectral index, T = 3.32. The 
upper limit, shown in Figure El is below the flux values 
from earlier in 2002, but higher than the average flux for 
the entire data set. 

The fit of the spectral index versus period to a constant 
yields a x 2 of 7.1 for 5 degrees of freedom, corresponding 
to a 0.21 x 2 probability, which is marginally consistent 
with being constant. Since the fit data points include only 
the statistical errors, the x 2 w iU be smaller with the in- 
clusion of systematic errors (^0.1 for the spectral index). 
Therefore, there is no evidence that the spectral index 
of PKS2155— 304 varies with time. However, the evidence 
also does not rule out possible temporal variability in the 



spectral index during the aforementioned large changes in 
flux. 

A correlation between harder spectra and 
higher flux values might be expected for PKS 
2155—304, since this trend was demonstrated for 
Mkn 421 |(Aharonian et al., 2002b)| and Mkn 50 1 
HDjannati-Atai et al., 1999| | Aharonian et al, 200 1| ). 
While Figure is suggestive of such a correlation, the 
inference is mainly drawn from the July 2002 data point, 
taken during the first month of H.E.S.S. operations 
where only one telescope was operational, and thus the 
fit spectrum is most likely to have larger systematic 
errors. In addition, the spectral index was shown to be 
marginally consistent with being constant. Therefore, 
further observations with the more sensitive complete 
H.E.S.S. array are needed to make a definitive statement 
regarding any correlation. 

Given the large significance observed in each month of 
data taking, a variability search on a monthly time scale 
represents a rather coarse temporal binning of the data. 
Therefore studies were performed on the data to determine 
if the VHE flux from PKS 2155-304 is also variable on 
shorter time scales. 



7.2. Nightly Time Scale 

The integral flux greater than 300 GeV measured from 
PKS 2155—304 is calculated for each night of observations 
and is shown in Figure [7] The exposure is different in the 
various nights and the spectral index assumed in the flux 
calculation is the value determined from the power law 
fit to the time-averaged spectra for 2003 (r = 3.32). This 
value is used since there is no strong evidence for variabil- 
ity in the spectral index. Table lists the results of the 
fit of the nightly fluxes during each month to a constant. 
The poor x 2 values for some of the dark periods show 
that variability in the flux of VHE gamma-rays from PKS 
2155—304 clearly occurs on time scales of days. However, 
the June and September 2003 results also show that there 
are periods of steady flux as well. 



Table 5. The x 2 , degrees of freedom (d.o.f.), and x 2 prob- 
ability, from the best fits to a constant of the nightly 
integral flux greater than 300 GeV observed from PKS 
2155—304 versus time for each dark period. 



Dark Period 


2 

X 


d.o.f. 


P(> x) 


07/2002 


6.4 


3 


9.5xl0 -2 


10/2002 


21.8 


8 


5.4xl0~ 3 


06/2003 


5.7 


7 


0.58 


07/2003 


30.1 


14 


7.4xl0~ 3 


08/2003 


74.9 


12 


3.7X10" 11 


09/2003 


1.0 


2 


0.60 
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Fig. 7. The nightly integral flux greater than 300 GeV versus time measured from PKS 2155—304 in a) July 2002, b) 
October 2002, c) June 2003, d) July 2003, e) August 2003, f) September 2003. Only the statistical errors are shown. 
The dashed lines represent the best x 2 fit of the data to a constant for each darkness period. 



7.3. Sub-hour Time Scale 

A search for evidence of intra-night variability in the VHE 
flux from PKS2155-304 is performed for the 2003 data 
set. A run by run (28 min) time scale is chosen for sim- 
plicity and because each stereo run should yield on average 
^4er given that ~44ct was observed in 55.4 live time hours 



of stereo observation. The study is not performed for the 
2002 data because the mono-configuration is considerably 
less sensitive and the data are more sparse (one on-source 
run on most nights). The respective fluxes greater than 
300 GeV for each run are fit to a constant for each night 
of observation, with the hypothesis that any variability 
will manifest itself in the form of a poor x 2 . Only runs 
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Fig. 8. a) The distribution of x 2 probabilities from best 
fits to a constant of the run by run integral flux greater 
than 300 GeV for each night of H.E.S.S. observations of 
PKS2155-304 in 2003. The solid line is for all nights of 
data taking and the dashed line represents nights with 
more than one degree of freedom (3 runs minimum), b) 
The run by run integral flux greater than 300 GeV ob- 
served from PKS 2155-304 for the night of MJD 52821. 
Only the statistical errors are shown. 



with zenith angles of less than 45 degrees are analyzed to 
eliminate possible systematic effects from large zenith an- 
gle observations. As motivated in the previous section, the 
assumed spectral index for the flux calculation is T — 3.32. 
Since there are differing numbers of runs in a given night, 
the x 2 probability is calculated from the fit to a constant 
for each night and is shown by the solid line of Figure|HIa) 
for all nights with a minimum of one degree of freedom. 
The dashed line in Figure |Sf a) represents the nights with 
more than one degree of freedom. A larger population of 
nights with low probabilities exists and is evidence for 
intra-night flux variability. 

Figure El[b) shows the flux light curve for the night 
with the lowest probability of being consistent with con- 
stant along with the corresponding fit (dashed line). On 
this night, MJD 52821, flaring behavior is suggested since 



the flux increases by a factor of 2.7±0.7 on the time scale 
of half an hour, and then decreases by a factor of 2.3±0.9 
in a similar time. Unfortunately at this stage in the con- 
struction of H.E.S.S., observations were only possible with 
the less sensitive two-telescope software stereo configura- 
tion, thus a more detailed study is not possible due to the 
relatively small number of observed gamma rays (~1 per 
min). However, with the completed H.E.S.S. experiment 
it will be possible to explore similar phenomena on a time 
scale of a few minutes. 

8. Conclusions 

Gamma rays with energies greater than 160 GeV have 
been detected with high significance (^45c) from the HBL 
object PKS 2155-304 in H.E.S.S. observations made in 
2002 and 2003. Although observations were made in less 
sensitive configurations (fewer telescopes) of H.E.S.S. be- 
fore the completion of the construction of the full experi- 
ment, a strong detection is found in each of the dark peri- 
ods of observations with appreciable exposure. The VHE 
flux observed from PKS 2155—304 is shown to be vari- 
able on the time scale of months, nights, and hours with 
an average flux above 300 GeV between 10% and 60% 
of that observed from the Crab Nebula. An extreme case 
of variability shows an increase in the observed flux by a 
factor of 2.7±0.7 in a 30 minute interval followed by a de- 
crease by a factor of 2.3±0.9 in the following 30 minutes. 
Energy spectra are made for each dark period for which 
a detection occurred and are found to be characterized 
by a steep power law. No evidence is found for variabil- 
ity in the spectral index of PKS 2155—304 over time, or 
for a hardening of the spectrum with increased flux lev- 
els. However, these hypotheses are also not ruled out. A 
time-averaged energy spectrum is determined for the 2003 
H.E.S.S. observations and fits to a power law with an ex- 
ponential cutoff (r=2.90t^23, E CVLt =lA±%% TeV), and a 
broken power law (ri=3.15±^, r 2 =3.79j;^, ^break = 
0.7 ± 0.2 TeV), show an improved \ 2 P er degree of free- 
dom relative to the one obtained from a fit to a simple 
power-law (r = 3.32 ± 0.05). Although the fits to power 
laws with features are not significantly better than that 
of a fit to a simple power law, the suggestion of features 
may indicate the absorption of TeV photons on the extra- 
galactic infrared background light. However, any spectral 
curvature could also be intrinsic to the blazar. 

Comparison of the properties of PKS 2155—304 mea- 
sured by H.E.S.S. to those of other blazars detected at 
TeV energies is somewhat difficult due to the highly 
variable nature of these objects and their differing red- 
shifts. Regardless, comparison of the flux (magnitude 
and variability) and spectra, may yield some insights 
into phenomena associated with TeV bright AGN. In 
the cases of most well studied TeV blazars, Mkn 421 
(z = 0.030) and Mkn 501 (z = 0.033), the observed flux 
has been found to vary from a fraction of the Crab flux 
to peak values of ~7 (Aharonian et al., 2002b) and ~10 
(Aharonian et al, 1999a) times the Crab flux respectively. 
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The flux detected from 1ES 19 59+650 (z = 0.047) ranges 
from 5% to 220% of the Crab [(Aharonian et al, 2003bT 



and the flux measured from 1ES 1426+428 has shown vari- 
ability by a factor of 2.5, with peak fluxes of order 10% 



of the Crab (Aharonian et al, 2003a) Clearly, the flux 
observed from PKS 2155—304, which ranges from 10% to 
60% of the Crab, does not vary as much as that from Mkn 
421, Mkn 501 and 1ES 1959+650. Additionally, the peak 
fluxes observed from these objects are considerably higher 
than that measured for PKS 2155—304. The observed 
flux from PKS 2155—304 is most similar, both in mag- 
nitude and variability, to that found from 1ES 1426+428. 
Comparing the spectra of the other AGN to that found 
from PKS 2155—304 also shows some differences. The 
spectral indices, prior to the exponential cutoffs found, 
of Mkn 501 (r = 1.92 ± 0.20), Mkn 421 (T = 2.19 ± 0.04), 
and 1ES 1959+650 (T = 1.83 ± 0.17) are harder than 
r = 3.32 + 0.05 found for PKS 2155-304. Again the most 
similar blazar is 1ES 1426+428, whose spectral index is 
T = 3.50 ± 0.35 |(Petry et al, 2002)| It is interesting that 

0.117) are most 
0.129) is 



the properties of PKS 2155-304 (z 
like those of 1ES 1426+428, whose redshift (z 
also the most similar. However, it cannot be stated from 
the evidence presented that the similarities are a direct 
consequence of the comparable redshifts. 

As of 2004, the H.E.S.S. detector is complete and is 
considerably more sensitive than any of the configurations 
presented here, allowing for more detailed studies of PKS 
2155—304 and other astrophysical objects to be performed 
in the future. 
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